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ABSTRACT 

The Hubble Deep Field South (HDF-S) Hubble Space Telescope (HST) observations are 
expected to begin in October 1998. We present a composite spectrum of the QSO in the HDF-S 
field covering UV/optical/near IR wavelengths, obtained by combining data from the Australian 
National University 2.3m Telescope with STIS on the Hubble Space Telescope (HST)[J This 
intermediate resolution spectrum covers the range 1600-10000A and allows us to derive some 
basic information on the intervening absorption systems which will be important in planning 
future higher resolution studies of this QSO. 

The QSO J2233-606 coordinates are 22:33:37.6 -60:33:29 (J2000), the magnitude is B=17.5, 
and its redshift is z em = 2.238, derived by simultaneously fitting several emission lines. The 
spectral index is a — —0.7 ±0.1, measured between the Lyman-a and Mg II emission lines. Many 
absorption systems are present, including systems with metal lines redwards of the Lyman-a 
emission line at z a b s — 2.204, 1.942, 1.870, 1.787 and a few very strong Lyman-a features at 
z a b s = 2.077, 1.928, without similarly strong metal lines. There is a conspicuous Lyman limit 
(LL) absorption system that is most likely associated with the z a bs — 1-942 system with a 
neutral hydrogen column density of ~Nhi = (3-1 ± 1.0) x 10 17 cm~ 2 . There is some evidence for 
the presence of a second LL absorber just to blue of the conspicuous system at z = 1.870. We 
have employed a new technique, based on an analysis of the shape of the observed spectrum in 
the region of the LL absorption, to explore the properties of the gas. We tentatively conclude 
that this system might have suitable characteristics for measuring the deuterium-to-hydrogen 
(D/H) ratio. 

Subject headings: quasars: absorption lines — quasars: emission lines 

1. Introduction 

Compared to the original HDF, scientific studies of the the Hubble Deep Field South (HDF-S) may 
benefit from the additional information provided by a bright high redshift z > 2 QSO lying close to (but 
not within) the targeted WFPC-2 imaging area. Initially 12 small (ldeg x ldeg) fields were selected as 
candidate target fields. Whilst no known high z QSOs were present in any of the possible fields, Hewett 
and Irwin found a high z QSO candidate in one of the fields using a UK Schmidt Telescope IllaJ objective 
prism plate scanned by the Automated Plate Measuring facility in Cambridge. Follow-up observations of 
the candidate at the Anglo-Australian Telescope confirmed the object to be a QSO (J2233— 606 22:33:37.6 
-60:33:29 (J2000), B=17.5) with a redshift z > 2 (see Boyle 1997 for details). We note for completeness 



1 Our ground-based 2.3m spectrum a nd the composite spectrum are available in electronic format at the address 
http: / /bat. phys.unsw.edu.au/~kms/hdfs/ 
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that 20cm radio observations have been made (Ray Norris, private communication reported in Boyle 1997) 
that give an upper limit on the radio flux as S20cm < 3mJy. There is a strong radio source 10 arcmin from 
the QSO (S 20cm = 146mJy). 

We present here a low-resolution spectrum of the HDF-S QSO taken with the Australian National 
University (ANU) 2.3m Telescope and combined with the publicly available STIS spectrum providing wide 
spectral coverage. The QSO has an optically thin Lyman limit which is a possible candidate for measuring 
the deuterium-to- hydrogen (D/H) ratio. The QSO also has several strong absorption features redwards of 
the Lyman-a emission line, to which we assign tentative identifications. 

2. Observations and Data Reduction 

The ANU 2.3m observation of the QSO was taken on 1997 October 21 at Siding Spring in clear 
conditions with s» 2 arcsec seeing. The Double Beam Spectrograph was used with a 300 line/mm grating in 
the blue arm (binned by 2 in dispersion), and a 158 line/mm grating in the red giving a resolution of 8 A in 
each arm. A dichroic was used to split the light at w 5400 A. Due to the difficulty in combining the red and 
blue spectra around the dichroic wavelength there is some uncertainty in our spectrum in the wavelength 
region A5400 - 5600A. The data were reduced using standard procedures with the IRAF ^ image analysis 
software. The resulting spectrum (A3300 — 10000A) is the weighted combination of 2 exposures (1200s and 
900s). 

We combined our ground-based 2.3m spectrum with the HST STIS spectrum taken on 1997 October 
30 with the G230L grating for 3700s (A1600 - 3100A resolution w 3.4A) and the G430L grating for 2200s 
(A2900 — 5700A, resolution «6A). Our 2.3m spectrum was not spectrophotometric (we used a narrow slit) 
so we scaled it before combining to have the same flux level as the STIS spectrum in the overlapping region 
of the two spectra. In order to preserve the higher resolution of the STIS data and because of uncertainties 
in the flux calibration of the ground-based data in the blue, we only used data redward of the Lyman-a 
emission (4170A) from the 2.3m data in the combined spectrum. A comparison of the 2.3m spectrum and 
the STIS spectrum in the overlap region made by cross-correlating the two spectra does not show any 
significant wavelength shifts, with a typical uncertainty of 40 km s _1 (0.5A). This is consistent with the 
rms residuals obtained for the comparison lamp exposures for the 2.3m blue DBS spectrum, which were 
around 0.6 A. 

The final composite spectrum (shown in Fig. 1) has been rebinned to the dispersion of the UV 
STIS component of the spectrum (A < 2900), 1.55A/pix. The signal-to-noise per pixel is ~ 5 below the 
Lyman-limit break (A < 2700), w 20 in the Lyman forest, and up to « 40 around A7000 - 8000A. 

3. Results 



2 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. 
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3.1. Emission lines and QSO continuum 

We measured the redshift of each of the emission lines in the QSO independently and also determined 
the QSO redshift by simultaneously fitting some of the emission lines using VPFIT (Webb 1987, Carswcll 
et al. 1992). Table 1 shows, for each of the emission lines, the redshifts, the equivalent widths (EW) and 
the velocity widths. The four lines used in the simultaneous measurement were C III] + Al III, C IV and 
Mg II. The individual redshifts for each line are listed in Table 1. The redshift from the simultaneous fit is 
z=2.238. The difference of 1390 km s" 1 between the (C IV and C III] + Al III) and Mg II redshifts is not 
unusual (Espey et al. 1989). The observed EWs lie in the ranges given by the composite QSO spectra of 
Francis et al. (1991) and Boyle (1990). 

The spectral index of the combined spectrum is a = —0.7 ± 0.1 (where f v oc z/*), consistent with 
values derived for optically selected QSOs (eg. the LBQS survey yields 0.0 < a < —0.8; Francis et al. 
1996). The measurement was made by fitting to three sections of the spectrum between A4000-9000A free 
of detectable emission lines. Our estimate of the uncertainty in a was derived by fitting the combined 
spectrum and the 2.3m data alone and hence is approximate. 

3.2. Absorption systems 

An initial investigation of the QSO reveals that there is complex absorption present from intervening 
systems. The absorption system details are tabulated in Table 2. There are four systems with metal lines 
and two systems with no detectable metals (at the signal-to-noise and resolution of this data) rcdwards of 
the Ly-a emission line. 

There is a Ly-a line at z= 1.870 (L3) with associated metal lines (CII, Silll and strong CIV). If NV is 
present, it is blended with the stronger (z=1.928) of the 2 Ly-a lines at z=1.928 (LI) and z=1.942 (L2) 
(see top panel of Figure 2) which we are associating with the LL break (discussed later). We note that 
the CIV absorption line associated with L3 (z= 1.870) line falls blueward of the SilV emission line so a 
SilV identification is also feasible. However, this seems unlikely because there is no corresponding Ly-a 
absorption line for a SilV identification. 

Identifications of other metal systems are less secure due to there being fewer metal lines. Prominent 
Ly-a absorption features at z=1.942 (L2), z=1.787, z=2.204, all appear to have CIV absorption. Also, 
there are absorption lines lying mid-point between the Ly-a and NV emission lines which may be due to 
NV absorption corresponding to the z=2.204 system. 

There appear to be no detectable metal lines associated with the strong Ly-a line (z=2.077), or with 
LI (z=1.928). The strong line at 3003. 9A is due at least in part to Ly-/3 associated with LI (z=1.928). 

If the identifications above are correct, we are seeing a rich complex of CIV absorption spread over the 
redshift range 1.787 to 2.204 corresponding to a scale ~ 100 Mpc, if interpreted as a spatial separation. 
Detailed studies of this QSO and the region of the sky around it, particularly photometric redshifts analyses, 
should reveal whether we have discovered a huge structure formed at an early epoch, or whether the line of 
sight to this QSO happens to intersect multiple isolated galaxies or clusters. 

The STIS spectrum can be used to explore the properties of the Lyman limit absorption system to 
determine the neutral hydrogen column density, N#/, and to constrain the velocity width of the absorbing 
gas. This is of particular interest from the point of view of seeing whether the system has suitable 
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characteristics for measuring the deuterium to hydrogen ratio (D/H). The important requirement, apart 
from having a sufficiently high Njji for Dl to be detectable, is that the DI feature is not blended with an 
ill-placed HI line, blue-shifted with respect to the main system. Since the DI feature is blue-shifted with 
respect to HI by —81 km s _1 , any absorption complex with velocity structure broader than this may not 
be useful for a D/H determination (unless by chance a high neutral hydrogen column density component 
happens to lie on the blue side of the complex, relatively well separated from other components) . To explore 
the extent of the velocity spread in the LL absorber here, we introduce a method that utilizes the shape of 
the spectrum at the Lyman limit to estimate the width of the complex. Lyman limit absorption systems 
which have a low b-parameter and a single high column density component will produce a sharp Lyman 
limit discontinuity. Conversely, a high b-parameter or complex absorption will produce a more gradual 
decline at the Lyman limit. 

We first estimate the neutral hydrogen column density (N#j) of the absorption system causing the 
depression at the Lyman limit using N HI a LL = t where a LL = 6.3 x 10~ 18 cto 2 and —lnl = t (I is the 
ratio of the flux in the spectrum to the flux in the continuum just bluewards of the drop in the LL (2683A) 
and r is the optical depth) . To estimate the approximate uncertainty in Njji due to errors in the placement 
on the continuum, we derived N#j for two different continua (the lower of these, the preferred one, is shown 
in Fig. 1). The value derived is N HI = (3.5 ± 0.3) x 10 17 cto~ 2 . 

There are additional uncertainties involved with our estimate of Njyj; a general background of Lyman-a 
and higher order absorption lines will also contribute to the flux decrement we measure below 2683A. At 
z ~ 2, this effect means the estimate of N#j above is systematically high by approximately 10-20%. 

To identify the Lyman-a line associated with the Lyman limit (LL) we used the wavelength (2685A) of 
the LL to estimate where the associated Ly-a line would fall. Fig. 2 shows there are two candidate lines, 
one at z=1.928 (LI) and the other at z=1.942 (L2). They are sufficiently close to each other that it is 
unclear which is associated with the LL. 

In order to explore whether LI or L2 is associated with the LL, we generated a series of simulated 
spectra at the column density derived above, for a range of velocity dispersion parameters, (b = 10, 30 
and 60 km s _1 , where b = y / (2)cr) at the two redshifts of LI and L2. We then compared the agreement 
between the model and observed spectra in the LL region. The low spectral resolution makes it difficult 
to derive useful information by comparing model and observed individual Lyman series lines. Instead, a 
more promising method is to look explicitly at the wavelength and shape of the LL drop, at and around the 
point at which the flux falls to roughly a constant level. The second and third panels of Figure 2 show the 
2 possible systems associated with the LL with the models overlaid. 

Fig. 2 shows that it is still difficult to unambiguously determine whether the actual Lyman limit 
absorber is associated with LI or L2. However, if it is associated with L2, the b=30 and 60 km s _1 models 
are inconsistent with the data and we can only get a reasonable fit to the LL for b«10 km s _1 . The LI 
system requires a higher effective b-parameter, suggesting that HI has complex velocity structure, so that 
it would not be suitable for a D/H study. 

A further possibility we considered was that neither LI nor L2 were responsible for the LL, but 
that L3, (see Figure 2), was responsible for the depression. In order for its LL to cut in at the observed 
wavelength, N#j would have to be sufficiently high such that the high order Lyman lines were strong 
enough to blend together forming a substantial drop in the transmitted flux redwards of 912A in the rest 
frame. However, this would require a high neutral hydrogen column density which would leave no residual 
flux at wavelengths shorter than the LL, which in turn would require the zero level of the spectrum to be 
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incorrect. We nevertheless explored this possibility further by using the EW of L3 to estimate and 
then comparing its associated LL with the observed spectrum. This procedure should provide an upper 
limit on Njji, since in reality L3 is probably composed of multiple components. We found that the 
derived in this way did not produce a large enough flux decrement in the wavelength region shortwards of 
the observed drop at 2680A. 

Careful examination of Fig. 2 also suggests that the residual flux in the data between 2650-2680A is 
slightly higher than the synthetic LL from L2 (or LI) alone. The spectrum appears to drop to a lower 
average flux below ~ 2640A implying the possible presence of an additional high cloud. 

This suggests two steps in the LL, one associated with L2 (z=1.942) and one with L3 (z=1.870). We 
again used VPFIT to model the data with a 2-component LL system, fixing the two redshifts (z=1.942 and 
z=1.870), allowing ~Nhi to vary. The resulting fit can be seen in panel 4 of Fig. 2. This procedure results 
in a good fit to the data, suggesting that the observed LL is associated with L2 and that there is indeed a 
further weak LL absorption associated with L3. 

The VPFIT estimate of the L2 column density is N HI = (3.1 ± 1.0) x 10 17 cm~ 2 . A 2-cr upper limit 
was also obtained for the N#j of system L3. This was determined to be Nhi = 3.2 x 10 17 cm~ 2 . If the LL 
is associated with L2 this system may turn out to be suitable for measuring D/H, because the (tentatively) 
low effective b-parameter suggests simple velocity structure and a small velocity dispersion. 

Using the two QSO continua described above, we estimated the parameter representing the flux 
decrement between the Lyman- (3 and the Lyman-a emission lines, = 1 — {fobs/ fcont) = 0.14 ± 0.04. 

4. Discussion 

We have reported the results of a study of a low/intermediate resolution spectrum covering the 
UV/optical/ncar IR of the HDF-S QSO J2233-606, formed by combining HST/STIS observations with 
higher wavelength data from the ANU 2.3m telescope. The intrinsic QSO spectral characteristics have been 
measured from this data. 

A study of the absorption properties have revealed multiple CIV absorption systems, suggesting that 
either a large structure or multiple galaxies or clusters intersect the sightline to this QSO. Structure in the 
LL absorption is probably caused by the presence of at least 2 high N(ffl) systems, one of which may have 
properties which make it suitable for studying D/H. 

Despite the low spectral resolution, we have been able to derive reasonably detailed characteristics of 
the absorption systems. Future high resolution spectra will of course yield far more information. However, 
it is already clear from our results that this QSO and the HDF-S studies of its environs will reveal a host of 
fascinating new information about the z~ 2 Universe. 
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Fig. 1— Composite spectrum of the HDF-S QSO based on ground-based ANU 2.3m data (4170 - 10000A) 
and the HST STIS spectrum (1600 - 5700 A). Also shown is one (the lower) of the two continua used in the 
initial determination of Njjj . 

Fig. 2. — Models of the Lyman limit absorption: the top panel illustrates the three candidate Lyman-a 
lines which may be associated with the LL. The 2nd and 3rd panels show the LL region with our models 
overplotted to determine the velocity structure of the system associated with the LL. N HI was calculated 
as described in the text, and 3 values of the b-parameter were used (A=10, B=30 and C=60 km s^ 1 ) for 
each of the 2 redshifts. From the fits we can see that either LI (z=1.928) with a high b-parameter or L2 
(z=1.942) with a low b-parameter will fit the data. The bottom panel shows the LL region with a model 
overplotted including the combined effect at the LL from systems L2 and L3. 
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Line z EW ofcs (A) FWHM (km s" 1 ) 

C IV 2.236 97 10370 

Al III + CIII] 2.238 104 12560 

Mgll 2.252 192 9900 

all 3 2.238 



Table 1: Emission line parameters measured using the composite ANU 2.3m and HST/STIS spectrum. 
The columns give the emission line identification, the redshift for each line, the observed equivalent width 
(EW ofcs ) and the velocity width (FWHM). 



A (A) 


a(X) 


ID 


^■abs 


EW obs (A) 


o-(EW) 


3388.5 


1.5 


Ly-a 


1.787 


1.7 


0.2 


4316.8 


1.5 


CIV 


1.787 


1.6 


0.2 


3463.8 


1.7 


Silll 


1.870 


2.5 


0.3 


3488.5 


3.0 


Ly-a 


1.870 


6.8 


0.5 


3830.1 


1.3 


CII 


1.870 


1.1 


0.2 


4444.2 


1.9 


CIV 


1.870 


4.0 


0.4 


3003.9 


3.2 


Ly-0 


1.928 


7.3 


0.5 


3559.2 


2.0 


Ly-a 


1.928 


5.1 


0.5 


3577.3 


1.6 


Ly-a 


1.942 


2.2 


0.3 


4558.7 


1.8 


CIV 


1.942 


3.3 


0.4 


3741.0 


1.9 


Ly-a 


2.077 


4.9 


0.4 


3895.2 


1.6 


Ly-a 


2.204 


2.0 


0.3 


3973.3 


0.9 


NV 


2.204 


0.9 


0.2 


4965.3 


1.7 


CIV 


2.204 


2.7 


0.4 



Table 2: Absorption line identifications and parameters, with 1-cr errors, for the redshift systems discussed 
in the text. EW b s is the observed equivalent width. 
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